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Graphene Superlattices (GSs) have attracted a lot of attention due to its peculiar 
properties as well as its possible technological implications. Among these characteris¬ 
tics we can mention: the extra Dirac points in the dispersion relation and the highly 
anisotropic propagation of the charge carriers. However, despite the intense research 
that is carried out in GSs, so far there is no report about the angular dependence of 
the Transmission Gap (TG) in GSs. Here, we report the dependence of TG as a func¬ 
tion of the angle of the incident Dirac electrons in a rather simple Electrostatic GS 
(EGS). Our results show that the angular dependence of the TG is intricate, since for 
moderated angles the dependence is parabolic, while for large angles an exponential 
dependence is registered. We also find that the TG can be modulated from meV to 
eV, by changing the structural parameters of the GS. These characteristics open the 
possibility for an angle-dependent bandgap engineering in graphene. 

PACS numbers: 73.21.Fg, 78.67.De 
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Bandgap engineering or band structure engineering is a term coined in the late eighties 
to refer to a powerful technique for the design of new semiconductor materials and devices.- 
This technique is based on the ability to modify the energy bands arbitrarily and to tailor 
them for a specific application. Among the typical tools to achieve bandgap engineering we 
can mention: doping, materials with variable gap, band discontinuities and superlattices. 
These tools can be used alone or in combination to obtain a particular band structure. For 
instance, superlattices are artificial periodic structures that can be created by alternating 
semiconductors with disimilar bandgaps, see Fig. [TJr. The difference of bandgaps of the 
constituent materials as well as the super-periodicity of the structure turn out in a periodic 
band-edge profile of barriers and wells, giving rise to allowed and forbidden energy bands, 
commonly known as energy minibands and gaps respectively, Fig. |TJo. Semiconductor su¬ 
perlattices present a plenty of physical effects that can be exploited technologically, among 
them we can find excitonic effects, miniband transport, WannierStark localization, Bloch 
oscillations, electric field domains, resonant and sequential tunnelling.— Specifically, super¬ 
lattices can be used as injection medium or as active region in the well-known Quantum 
Cascade Laser.- 

Within this context, bandgap engineering in graphene is not the exception, and from the 
very beginning of the discovery of graphene^ - - the scientific community has tried to figure 
out how to create and modify a bandgap in graphene.- - — This in part due to the gapless 
dispersion relation of graphene^— as well as a natural technique to tailor the band structure 
of graphene-based structures for a specific application. Among the different approaches to 
create a bandgap we can mention graphene nanoribbons,- - - epitaxial substrates,— - — bi¬ 
ased bilayer graphene,- 3,11 hydrogen absorption,— 1 ^ graphene nanomesh,—— and strain 
engineering.— 1 ^ All these approaches are based on the aperture and modulation of a 
bandgap. ffowever, as we above mention there are another possibilities, for example su¬ 
perlattices. Indeed, superlattices in graphene have been studied intensely in the past few 
years.— - — Among the most important characteristic of graphene superlattices we can found: 
additional Dirac cones in the energy-dispersion relation and highly anisotropic propagation 
of charge carriers. The latter is quite interesting, since depending on the transversal wave 
vector or the angle of the electrons that impinging on the superlattice structure the propa¬ 
gation properties can be tuned readily. This opens the way for possible electron wave filters 
and laser devices, ffowever, as a first step to propose this kind of devices, it is important 
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to know in detail the angular dependence of the energy minibands and gaps, aspect that as 
far as we know is not reported up to date. 

Here, we propose a bandgap engineering based on the angular dependence of the prop¬ 
agation of Dirac electrons in graphene superlattices. Our approach relies on the formation 
of transmission minibands and gaps at diferente energy scales depending on the angle of 
incidence of Dirac electrons. In particular, we show that the angular dependence of the 
TG is intricate, being parabolic for moderate angles and exponential for large ones. To 
this respect, TGs from meV to various eV can be obtained by angular selection as well as 
by tailoring the structural parameters of GSs (number of periods, height of barriers, and 
widths of barriers and wells). So, angle-dependent bandgap engineering gives the possibility 
to modify the energy minibands and gaps in GSs almost arbitrarily and to tailor them for 
a specific application. 

The system we are interested in is a rather simple electrostatic graphene superlattice. 
In Fig. []J we show a schematic representation of this structure, it consists of a graphene 
sheet sitting on a non-interacting substrate like Si 02 , a back gate, and top gates arranged 
periodically along the superlattice axis (x). In the left part of the structure, it is also shown 
an electron impinging at a certain angle with respect to the superlattices axis. We have 
chosen this structure because from the experimental standpoint is more reliable. 46 ' 4 ' For 
instance, the energy of the incident electrons and the width and height of the potential 
barriers can be controlled by the doping of the substrate and the strength of the voltages 
applied to the back and top gates. The band-edge profile of this structure can be modelled 
as a periodic arrangement of potential barriers and wells, that in turn results in energy 
minibands and gaps, which strongly depend on the angle of the impinging electrons, see Fig. 
mi. The transmission properties of this system can be computed straightforwardly using the 
transfer matrix approach. 48,19 The basic information needed to apply this methodology is 
the dispersion relation, wave vectors and wave functions in the barrier and well regions as 
well as in the semi-infinite left and right regions. 50,51 In the well and semi-infinite regions 
the dispersion relation and wave functions comes as: 

E = ±Twpk, 


and 
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^±ik x x+ik y y 


( 1 ) 


where Vp is the Fermi velocity, k is the magnitud of the wave vector in these regions, k x 
and k y are the longitudinal and transversal components of k, and u± = sign(E)e ±ie the 
coefficients of the wave functions that depend on the angle of the impinging electrons, 
9 = arctan (k y /k x ). In the barrier regions these quantities come as: 


E = Vq ± Uvpq, 


and 



where Vo is the strength of the electrostatic potential, q is the magnitud of the wave vector 
in the barrier regions, q x and q y are the components of q, and v± the coefficients of the wave 
functions. 50,51 Once these quantities are known, we can apply the continuity conditions of 
the wave function along the superlattice axis as well as the consevation of the transver¬ 
sal momentum= q y ), and define the transmission probability in terms of the so-called 
transfer matrix, 



which depends on the transfer matrices of barriers and wells, and the number of periods as 


well, for more details see Refs. 50 and 51. 

At first, we want to discuss the formation of energy minibands and gaps as a function 
of the angle of incidence. In Fig. [2] we show the transmission probability as a function 
of the energy for three different angles. Specifically, Fig. [2^ depicts the transmittance for 
normal incidence and 9 = 5°, solid-red and solid-black lines respectively, while Fig. [2)3 
shows the case of 9 = 15°. The number of periods, the widths of barriers and wells, and 
the height of the barriers are N = 10, ds = 10a, dw = 10a and Vq = 1.0 eV, respectively. 
These parameters will remain fixed, virtually, throughout the study. As we can see Klein 
tunneling (perfect transmission) prevents the formation of energy minibands and gaps for 
normal incidence, 9 = 0°. Once the angle of incidence is different from zero, transmission 
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minibands and gaps start to develop. For small angles, as in the case of 5°, we have pseudo 
minibands and gaps,— since them are not well defined yet. By increasing systematically 6 
we will find that the mentioned pseudo minibands and gaps become well defined ones. This 
process takes place in staircase fashion, due to not all transmission minibands and gaps 
become well defined at the same angle. For instance, in the case of 9 = 15° we can see that 
the first and second gaps are almost well defined as well as the second and third minibands, 
however what we labeled as first miniband in reality is a pseudo miniband, since in the low 
energy side of it there is not a well defined gap. Even more, this miniband is shared among 
electrons and holes, positive and negative energies, respectively. By increasing the angle 
of incidence to 30° and 50° we can see that the widths of the energy minibands diminish, 
while the corresponding ones to the transmission gaps increase, see Fig. and b. In the 
specific case of the minibands, we can also notice that the number of resonances within them 
diminish as well. This reduction is quite important, since for example a particular miniband 
will occlude in a specific angle, and consequently the corresponding transmission gap will 
increase substantially. For instance, in Fig. [3b we can see that the first transmission gap is 
huge as a consequence of the collapsing of various minibands, or in other words the second 
miniband is at quite different energy scale, due to the occlusion of various transmission 
minibands. We can have a better perspective if we taking into account the contour plot 
of the transmittance as a function of both the energy and angle of the incident electrons, 
see Fig. [4]a. In fact, the contour tells us that the energy minibands and gaps are formed 
in a staircase fashion, the energy minibands have a semi-circular form, akin to a whisker, 
and are wider for small angles and virtually disappear for large angles. In addition, as 
the energy increases these whiskers bend, and the bending is steeper for higher minibands, 
as a consequence higher minibands occlude at lesser angles than lower minibands. Other 
important characteristics that we can find in EGSs are: 1) the number of resonances within 
a miniband is proportional to the number of periods in the superlattice, 2) by increasing the 
number of periods it is also possible to obtain well defined minibands and gaps irrespective of 
the angle of incidence,— except normal incidence, 3) by changing the widths of the barriers 
and wells as well as the height of the barriers it is possible to tune the number of minibands 
and gaps, and their energy location, see Fig. [4jo. 

Within this context, a natural question arises, what is the angular dependence of the 
transmission minibands and gaps?. In order to answer this question, and without lose of 
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generality, we will focus on what we labeled as first transmission gap. As we can see from the 
transmission contours, it is possible to infer that there are two angular regions for the first 
transmission gap due to the dramatic change that is presented about 60°-70°, see Fig. ED 
Indeed, we find a parabolic dependence in the range of 20° to 65°, while in the interval of 65° 
to 90° an exponential dependence is presented, see Fig. 0 Excellent parabolic and exponen¬ 
tial fits were found, see solid-red lines in Fig. [5j For instance, the Adj. R-Square is 0.99418 
and 0.99507 for the parabolic and exponential fit, which means that the fitting is quite good. 
The fitting parameters in the parabolic region are a = 6.93 x 10 -4 eV, b = —0.026 eV and 
c = 0.594 eV, while in the exponential region we have A = 2.99 x 1CD 10 eV, B = 0.322 0-1 
and C = 16.09 eV. It is also important to highlight that choosing appropriately the angle 
of incidence, the first transmission gap can be tune from meV (small angles) to various 
eV (large angles). This opens the possibility of a bandgap engineering with a tremendous 
energy range of modulation, notice the vertical scale in Fig. [5j This is a unique property 
presented by EGSs, since to modulate the bandgap in traditional semiconductor superlat¬ 
tices, it is mandatory to change mainly the constituent materials. Here, on the contrary, it 
is required to have full control of the angle of incidence, or more realistically speaking, total 
control over a narrow angular interval. To this respect, we can expect that with the current 
experimental technologies this requirement be totally affordable. In fact, recently, impor¬ 
tant steps have been made in order to find out the angular contribution of Dirac electrons 
to the transport properties in graphene-based devices.—Specifically, some groups have 
been able to discriminate the angular contribution by using tilting metallic electrodes (top 
gates) in single-barrier graphene structures. Finally, it is important to remark that in order 
to angle-dependent bandgap engineering be realistic and reliable other relevant issues have 
to be addressed. For instance, in the case of laser devices, studies about tunneling rates, 
relaxation times, tunneling escape probabilities, emission lifetimes and radiative efficiency 
are needed in order to elucidate the conditions for population inversion. 

In summary, we have proposed an angle-dependent bandgap engineering in graphene. 
This approach is based on the angular dependence of the propagation of Dirac electrons 
in EGSs. Specifically, transmission minibands and gaps can be formed at different energy 
ranges depending on the angle of the incident electrons. In the case of the transmission gaps, 
the main gap shows a parabolic and exponential angular dependence for moderate and large 
incident angles, respectively. The transmission gap can be tuned from meV to eV by angular 
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selection as well as by changing the structural parameters of GSs. So, in principle GSs offer 
a tremendous possibility to design and develop devices based on angle-dependent bandgap 
engineering. 
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FIG. 1. Schematic representation of Typical Semiconductor and Electrostatic Graphene Superlat¬ 
tices. The cross-sections of these superlattices are shown in (a) and (c), while in (b) and (d) the 
conduction band-edge profiles and the energy minibands (Ai and A 2 ) of them are depicted. cIb, 
d\y and Vo represent the widths of barriers and wells, and the height of the barriers, respectively. 
Aside to the structural differences, the main physical difference is that energy minibands and gaps 
in the case of Electrostatic Graphene Superlattices depend on the angle of incidence. 
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FIG. 2. Formation of energy minibands and gaps in EGSs as a function of the angle ( 9 ) of the 
impinging electrons. In (a) the transmittance for 0 = 0° and 9 = 5° is presented, while in (b) the 
transmission probability for 6 = 15° is shown. We have denoted the energy minibands as Ai, A 2 
and A 3 , and the energy gaps as AEgi and AEg 2 . 
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FIG. 3. Evolution of the energy minibands and gaps in EGSs for different angles of incidence: (a) 
30°, (b) 50° and (c) 80°, respectively. In general, the energy minibands become narrower, while 
the energy gaps get larger as the angle of incidence increases. Even more, a dramatic change is 
presented in (c), where the first transmission gap changes hugely. 
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FIG. 4. Contour plots of the transmittance as a function of the energy and angle of the incident 
electrons: (a) Transmission contour of EGS for superlattice parameters: N = 10, ds = dw = 10a 
and Vo = 1 eV, (b) The same as in (a), but with double of the size of the width of the quantum 
wells, dw = 20a. White arrows help us to guide the evolution of the first transmission gap, and 
specifically to indicate that at certain angle there is a dramatic change due to the occlusion of 
various transmission minibands. 13 
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FIG. 5. Angular dependence of the first transmission gap AEgi in EGSs. (a) In the angular range 
of 20° to 65° a parabolic dependence is presented, while (b) an exponential dependence arises 
from 65° to 90° . In particular, the parabolic and exponential dependence give the possibility to 
modulate AEgi from meV to eV. The solid-red lines correspond to the parabolic and exponential 
fit, respectively. The Adj. R-Square is 0.99418 and 0.99507 for the parabolic and exponential fit, 
indicating that the fitting is quite good in both cases. 
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